oo 



< 



o 



Semiclassical limit for generalized KdV equations 
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CN ' Abstract 

We study the semiclassical limit of the (generalised) KdV equation, 

for initial data with Sobolev regularity, before the time of the gradient 

catastrophe of the limit conservation law. In particular, we show that 

p^ I in the semiclassical limit the solution of the KdV equation: i) converges 

in if ^ to the solution of the Hopf equation, provided the initial data 
belongs to H", ii) admits an asymptotic expansion in powers of the 
semiclassical parameter, if the initial data belongs to the Schwartz 
class. The result is also generalized to KdV equations with higher 



rS^ • order linearities. 

-4— » 

Introduction 

psj , We consider the following class of partial differential equations 

>. 

^' dtu = a{u)dxU + ^eid'^''^'^u, u{x,0) = (p{x), (G) 

^' 

^—i . depending on a family of parameters e = (ei, . . . , £„) G R". Here a{u) 

t^ ' is a smooth function, a: G R, and the initial data is independent of e. 

We call the above family of equations generalized KdV equations, for 
it contains as a particular example the KdV equation itself: 

dtu — u dxU + e d^u . 

^^ ' Other examples are given by the Kawahara equation |15| . which is 

M , obtained by choosing a{u) ~ u and n = 2, as well as nonlinear gener- 

alizations of KdV, for n — \ and a arbitrary. 

We are interested in the behaviour of the solutions of ([G]) as the 
parameters e vary. In particular, we consider the behaviour when 
e -^ 0, in which limit, formally, ([G)) becomes a quasilinear conservation 
law. of the form: 

dtu = a{u)dxU, u{x,0) ^ ip{x). (H) 
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The case a{u) = w is known as Hopf equation. The problem of study- 
ing solutions of equation ([G|) as £ — >■ is known in the literature as the 
semiclassical (or singular, or dispersionless) limit. Since solutions of 
equation plj) may develop a singularity at a critical time 0^ < ic < oo, 
we study local-in-time solutions of (jG]) for those classes of initial data, 
to be specified below, for which the Cauchy problem for the generalized 
KdV equation is locally well-posed uniformly with respect to e G M". 

There is an extensive literature on the initial value problem for 
generalized KdV equations, which is based essentially on two distinct 
approaches: the first makes use of the inverse scattering and Riemann- 
Hilbert methods, and it applies to those equations of the class ([G]) 
which are integrable, such as KdV. The second applies to a more gen- 
eral class of equations, and makes use of fix-point arguments for the 
associated integral equation 

u{t]e) = W{t;£)u{0]e)+ W{t - s]e)a{u{s;e))u^{s;e)ds , (W) 

Jo 



W{t■,e)^exp{tJ2e^^'J^') 



Since the seminal paper of Lax and Levermore |18) . the method of 
inverse scattering has been successfully used to study the semiclassical 
limit of the KdV equation, both before and after the critical time ic 
of the Hopf equation. For time smaller than tc, rigorous results have 
been obtained for those initial data whose scattering transform can 
be computed in the semiclassical limit using the WKB analysis. In 
this case, the corresponding solutions are proven to converge in L? to 
the solutions of the Hopf equation [THl [22] . If, in addition, the initial 
data satisfy some analyticity assumptions, then the powerful nonlinear 
steepest-descent analysis [7| applies to the study of the semiclassical 
limit. For these initial data, the solutions are known to converge uni- 
formly - see 013]. 

The integral equation (W) has been a major topic of investigation 
since it was used by Kenig, Ponce and Vega [TC] to establish local 
well-posedness of ([G|) for polynomial nonlinearity a(u) and dispersion 
£ 7^ 0. In the particular case of the KdV equation, the authors of 
|16j have been able to prove local well-posedness in 7J*,s > |. Their 
results were then refined by many authors to obtain local and global 
well-posedness for low-regular initial data, see for instance O [T7J [S] . 

The theory of equation (W) relies heavily on the dispersive charac- 
ter of equation ([G]) , being based on the smoothing effects of the linear 
evolution operator W . Hence, it may be very difficult to apply it to the 
study of the semiclassical limit. Moreover, since the equation (H) may 
be ill-posed for s < |, it seems unreasonable to study the semiclassical 
limit for low- regular solutions. 

Due to the limitations of the inverse scattering transform and of 
the integral equation (W) , in order to deal with the semiclassical limit 



of the general equation ([G|) we choose a different approach, namely 
Kato's theory of quasi- linear equations |13| . This approach turns out 
to be particularly suitable to our problem because it allows to treat 
equation IG]) for all values of e G M" on the same footing, and it is 
very robust under perturbations|^ 

Following Kato, we consider ([G]) as a quasilinear equations on a 
Banach space X, of the form 

'^=A{u{e)-e)u + f{u{e)-e), 0<t<T, u(0) = ^. (Q) 

Here A{y; e) is a linear operator, depending on e e M" and on some 
element y £ X. In addition, for any fixed e and y the operator A{y] s) 
generates a Cq— semigroup on X . In our case, 

n 

Aiy;e)=a{y)d,+Y,£^€^\ / = 0, 

and we choose X = L^ (R) . Kato himself used his theory to construct 
local-in-time solutions of equation ([G]). In particular, he established 
local well-posedness for the KdV equation in H'^,s > |, both when 
£ ^ [14] and when e = [13]. However, he did not consider the 
semiclassical limit. 

We establish simple conditions, under which the local-in-time solu- 
tion of the Cauchy problem ((G)) with initial datum in H^ is continuous 
and A^— differentiable with respect to e. 

Essentially, we show that: 

• If s > 271 -I- 1, the local-in-time solution u{e) of (G) is continuous 
with respect to £ e M" . 

Let K = J2"^^ N,{2i + 1) and N = J2'Li ^^- U s-K >2n + l, 
then the partial derivative 



• 



9£f ...a£^ 

exists in H''^^ and it is continuous with respect to £ e 



The above results can be applied to the project, proposed by Du- 
brovin and Zhang |10l [8] , of Hamiltonian perturbations of quasilinear 
conservation laws (H) . Indeed, any equation the form (H) can be writ- 
ten as an infinite dimensional Hamiltonian system; within this theory, 
one looks for suitable deformations - depending on arbitrary functions 
of u and its derivatives - such that the equation remains Hamiltonian. 
One simple example is given by KdV, which can be obtained as a 
Hamiltonian perturbation of the Hopf equation. 



^An alternative method, which may be suitable for the study of the semiclassical limit, 
is given by the Bona-Smith energy method [2]. 



In addition to the above problem, the project includes a classifica- 
tion of integrable perturbations and a characterization of the solutions 
of the perturbed equations, both before the critical time and in a neigh- 
borhood of it. It should also be noted that the project applies not only 
to single equations, but also to systems of first oder quasilinear PDEs, 

Before critical time, the Dubrovin-Zhang theory provides a way to 
construct solutions of the perturbed equation in terms of solutions of 
the unperturbed one. Without going into detail, one looks for solutions 
of the perturbed equations as formal power series u = v'^ + v^ e'^ + . . . , 
and argues that i;*^ is a solution of the unperturbed equation, while 
the subsequent coefficients can by obtained from u" by a recursive pro- 
cedure. This construction, although extremely powerful in predicting 
behaviour of the solutions, is based on formal identities, and requires 
rigorous justification. 

The paper is organized as follows: after recalling in Section [T] some 
basic elements of the theory of Cq— semigroups, we consider in Section 
[5] the results obtained by Kato in [T3], which we will use in order to 
prove our results. 

Sections |3] and H] are the core of the paper. We first prove the exis- 
tence of a positive time T for which the solutions of the problem ((G)) . 
with initial data in H'^, s > 2n + 1, are continuous functions of the 
parameters e. In particular, such solutions are continuous as e — > 0, 
implying iJ'*— convergence to the solution of (H) in this limit. This 
result generalizes the one obtained by Lax and Levermore (before the 
critical time) to equations of type ((G|) - which are not necessarily in- 
tegrable - and to initial data in the Sobolev space H^. 

In Section H] we consider differentiability of solutions of ([G| with 
respect to e. Although our result - with suitable modifications - holds 
for every equation of type ([G|) , for simplicity we consider in detail the 
KdV example only. We show that if the initial datum of the KdV 
equation lies in the Sobolev space H^ , then the solution of the Cauchy 
problem is A^— times differentiable with respect to e, for N = [s/3 — Ij . 

In Section[5]we present the Dubrovin-Zhang theory of Hamiltonian 
perturbations of equation (H) , considering those aspects of the theory 
which are directly related to the results of the present paper: the clas- 
sification results of Hamiltonian perturbations and the construction of 
the solutions of the perturbed equations before critical time. 

In the last Section we apply the results obtained in Sections [3] and 
m to the Dubrovin's theory. Since equations of type ([G]) can be seen 
as Hamiltonian perturbations of (H), we provide - for this class of 
equations - a rigorous justification to the heuristic results of Section 
[5] In addition, we find an explicit formula for the coefficient v^ of the 



solution of a generic Hamiltonian perturbation of (H) in terms of the 
solution v'^ of the unperturbed equation. This is of the form: 



2dx\ (l + ia'i;0)2 

where a = a(v^) is the non-linearity of (H), and c = c{v^) is a function 
characterizing the Hamiltonian perturbation. 

Notation 

Given the real Banach spaces X,Y, . . . , we let [| ||^- , || [|y , . . . de- 
note the corresponding norms. C{Y, X) denotes the Banach space of 
bounded linear operators from Y to X, with norm || \\y x^ while C{X) 
denotes the Banach space of bounded linear operator from X to it- 
self with the norm || ||^. We call D{A) the domain of an operator 
A. L^ :— L^(R) denotes the Hilbert space of square integrable real 
functions and H'^ :~ iJ'*(K), s > denotes the Sobolev space of order 
,s. The symbol d" denotes the n-th derivative with respect to x or the 
corresponding operator on L^ with domain H". 

In the present paper we consider real Banach spaces and real func- 
tions only. 
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1 Cq— semigroups 

The theory of Cq— semigroup is a standard tool in analysis. Following 
[13) and [20], we recall the elements of the theory we will use in the 
rest of the paper. 

A one parameter family of linear operators {T{t), <t < oo} on 
a Banach space X is a Cq — semigroup if it is a strongly continuous 
semigroup of bounded linear operators, namely, it satisfies: 

• T(0) = /, T{t)T{s) = T{t + s), i,s>0, 

• limtj_o T{t)x = X, V.T G X . 

Here / is the identity operator on X. The linear operator defined by 

f T(t)x — X 

D(A) = <^ a; e X : lim —-- exists 

1 tio t 



and 

Ax = lim ^^^^fc^, Vx e D{A), 

is the infinitesimal generator of the Cq— semigroup {T(f)}. The oper- 
ator A is closed and densely defined. 

A standard theorem shows that for any Cq— semigroup there exist 
two positive constants M >1 and /3 > such that 

\\T{t)\\x<Me^\ 0<t<oo. (1) 

In particular, a Cq— semigroup with constants Af = l,/3 = is called 
a semigroup of contractions. We denote by G{X, M, f3) the set of in- 
finitesimal generators of Cq— semigroups with constants M, (3. 

Remark 1. Note that if A e G{X, l,/3), then A- [il e G{X, 1,0). 

Later on we will need a perturbation theorem for generators of 
semigroups of contractions; for this purpose we introduce the following 
notions. 

An operator ^ on a Hilbert space X is said to be dissipative if for 
every x £ D{A), we have Re{Ax, x) < 0. 

If A and B are operators on a Banach space X, we say that B 
is relatively bounded with respect to A with relative bound p > if 
D{A) C D{B) and there exists a cr > such that 

\\Bx\\^<p\\Ax\\^ + cy\\x\\^, yxeD{A). (2) 

Theorem 1. Let A G G(X, 1, 0) he the generator of a Ca^semigroup 
of contractions on a Hilbert space X . Let B be dissipative and relatively 
bounded with respect to A with relative bound p < \. Then A -\- B is 
the generator of a semigroup of contractions. 

Proof See gD], Corollary 3.3 D 

The definition of dissipative operator and Theorem [T] can be gener- 
alised to any Banach space, with slight modifications. However, such 
generalisation is not necessary in our study. The following examples 
will be useful in the rest of the paper. 

Example 1. Let X be a Banach space and A be an anti-self-adjoint 
operator. Due to the Stone Theorem, A generates a Co^group of uni- 
tary operators, hence a Cq — semigroup of contractions. For instance, 
the derivative operator d^"'~^^ , n S N, with domain Z?(9^"+^) — H^"'^^ , 
is an anti-self-adjoint operator on the space L^ . More generally, the 
operator 

n 

Z?f+i=^£,9|'+\ e= (£!,...,£") eM", (3) 

with domain ij2)i+i j^ anti-self-adjoint on L^ for any value of the pa- 
rameter e. 



Example 2. J13\/ Let f{x) be a bounded differentiable function with 
bounded derivative on the whole real axis, take X = i^(]R), B = f{x)dx 
and D{B) = H^ . The operator B can be decomposed into an anti-self- 
adjoint part and a bounded self-adjoint part, B = Bi + B^ with: 

Bi = (/O, + i/,), B2 = -i/. . 

In particular, Bi has domain H^ and is anti-self-adjoint, thus, it gen- 
erates a Cq — semigroup of contractions. On the other hand, B2 is a 
bounded self-adjoint operator with norm 115211^2 = 5 ^up^gjj |/j;(a;)|, 
and a simple computation shows that B2 — ||i?2|li2 / is dissipative. Due 
to Theorem]^ we have that B — |[i?2 11^,2 1 generates a Cq — semigroup 
of contractions. Therefore, B G G{X,l,/3) with /3 = ||i?2||^2- In ad- 
dition, B is relatively bounded with respect to 9^""*"^, n > 1 with any 
relative bound p > 0. Due to Theorem]^ and the above discussion, 

B + edl''+^ eG{X,l,p), 

for any e G R. 

Example 3. Let g{x) be a continuous bounded function on M, take 
X = L^(]R), and let C = g(x) be the operator of multiplication by 
g, with D{C) ~ L^ . We have that C is a bounded operator, with 
norm j3' = supx(^K\g{x)\. Due to Theorem]^ if A G G{X,1,P) then 

A + CeGiX,l,p + p'). 

2 Kato's theory on quasilinear equations 

In this section we review Kato's results on quasilinear equations on a 
Banach space X, of the form 

du 

— ^A{t,u)u + f{t,u), 0<t<T, u(0) = (^. (4) 

Here A{t] y) is a linear operator, depending on the time t and on some 
element y & X , and such that for any fixed t and y the operator Ait; y) 
generates a Cq— semigroup on X. It should be noted that we do not 
present these theorems in their strongest form, but in a form adequate 
to our purpose. For the reader's convenience, we follow - as much as 
possible - the notation of the original paper [T3] . 

First, we consider the linear case, when the operator A and the 
forcing term / do not depend on u. 

Theorem 2. The linear non-homogeneous Cauchy problem 

du 

— ^A{t)u + f{t), 0<t<T, 7.(0) er, (5) 

has a unique solution 

u{t) e C{[0,T];Y) nG\[0,T];X), 
provided the following assumptions are satisfied: 



(i) X is a Banach space and Y C X is another Banach space, con- 
tinuously and densely embedded in X. Moreover, there exists an 
isomorphism S of Y into X. 

(a) There exists a positive j3 such that A{t) € G{X, l,/3) Vi G [0, T]. 

(Hi) SA{t)S^^ - A{t) = B{t) G C{X), and t ^ B{t) is a continuous 
operator valued function. 

(iv) Y C D{A{t)) so that A{t) G C{Y,X), and t ^ A{t) is a continu- 
ous operator valued function. 

(v) u(0) eY and f e Ci[0,T];Y). 
Proof See [H], Theorem I and II. D 

The following is a perturbation theorem for the linear equation ([5]) . 

Theorem 3. In addition to the assumptions of Theorem\^ consider 
the sequence of Cauchy problems 

dv" 

— =A^{t)u^ + r{t), 0<t<T, w"(0)Gy, (6) 

and suppose that, for any fixed n, the operator A" satisfies conditions 
(i) through (v) of Theorem\^ Moreover, suppose that 

(vi) A^{t) -^ A{t) strongly in C{Y,X), and sup^^rQ^,! ||^"(Ollyx ^^ 
uniformly bounded in n. 

(vii) B"'{t) — >■ B{t) strongly in C{X), and siipffriQ x] \\B"'{t)\\-^ is uni- 
formly bounded in n. 

(viii) it"(0) -^ u(0) in Y and f"- ^ f m C{[0,T],Y). 

Then M"(t) -^ u{t) in C{[0,T],Y) n C^{[0,T],X), where u"(t) is the 
unique solution of (^ and u{t) is the unique solution of Q). 

Proof See [H] Theorem V-VI. □ 

We now move to the analogue results for quasilinear equations of 
the form ^. For our purposes, it is sufficient to consider only the 
homogeneous case, when / = 0, namely: 

du 

— = A(t, u)u, < i < T, u(0) eW CY . (7) 

dt 

Here the set M^ is a bounded subset of Y. Due to the nonlinearity, 
existence of the solutions is not guaranteed on the whole time interval 
[0, T], but ~ in general - only for a smaller time T', with <T' <T. 
The reason we restrict to a bounded subset W is because we expect 
the time of existence to depend on the norm of the initial data. 

Following Kato, we make the following assumptions: 

(X) X is a reflexive Banach space and Y C X is another reflexive 
Banach space, continuously and densely embedded in X. There 
is an isometric isomorphism S of Y into X. Moreover, we fix a 
ball W CY oi radius R and centered in 0. 



(Al) There exists a positive (3 such that A{t^y) G G(X, 1,/?), for all 
f e [0,T] and 2/ £ M^. 

(A2) For any t, y G [0, T] x W, we have 

5'A(i)S-i - A{t) = B(i) G /:(x), 

and \\B{t)\\x < Al. 
(A3) For all t,y G [0,T] x W^, we have A{t) e C{Y,X). Fixed y eW, 
the function f n> A(t, y) is a continuous operator-valued function, 
and fixed t G [0,r], y —^ A{t,y) is Lipschitz continuous, in the 
sense that there exists a /ii such that 



\\A{t,y)~Ait,z)\\y^^<^l,\\y-z\\ 



X 



Theorem 4. Suppose conditions (X),(A1),(A2), (A3) are satisfied. Then, 
the quasilinear homogeneous Cauchy problem (^ has a unique solution 

u{t) e C{[0,T'];W) nC\[0,T'];X), 

for some < T' < T. In addition, T' has a lower bound uniquely de- 
pending on /3,Xi, ^i,R and monotonically decreasing in each variable. 

Proof. See ^ Theorem 6. □ 

We now state the perturbation theorem in the case of quasilinear 
equations of type ([7]) . 

Theorem 5. In addition to the assumptions of Theorem[A consider 
the sequence of quasilinear homogeneous Cauchy problems 

dii"- 

-- = A"(i,u")u", 0<t<r, w"(0)gW', (8) 

dt 

and assume that conditions (X), (Al), (A2) and (A3) of TheoremVM 
are satisfied for every n, with constants /3, Ai,^i independent on n. 
Moreover, suppose that 

(A4) \\B'^{t,y) - B''{t,z)\\^ < ^i2\\y - z\\y, uniformly mn. 

(CI) ^"(i) -^ A{t) strongly in C{Y,X). 

(C2) B'^it) -^ B{t) strongly m C{X). 

(C3) u"(0) G W and m"(0) -^ u{Q) in Y , as n^ oo. 

Then, there exists a positive time < T" < T , such that there is a 
unique solution 

u" G C([0, T"],W) n C^ ([0, T"],X) 

of 0), for every n, and a unique solution u of ^ in the same class. 
Moreover 

u"{t)^u{t) in C{[0,T"],W)nC^{[0,T"],X), 

as n i—> oo. The time T" has a lower bound uniquely depending on 
/3, Xi, fii, fi2, R, and monotonically decreasing in each variable. 

Proof See [B], Theorem 7. □ 



3 Continuity of solutions for generalised KdV 
equations 

Here we apply Kato's theory to study local-in-time solutions of the 
following family of Cauchy problems: 

nil 

— =A{u;e)u, u{0)&H^, (9) 

n 

A{u; e) = a{u)d,. +Y,£^€^' ^ 

depending on a family of parameters £ = (ei, . . . , e„) e M". Here a{u) 
is a smooth function, the initial value u(0) is independent on the Si for 
all i, and s > 271 + 1. We call the above equations generalized KdV 
equations. We are interested in the behaviour of the solutions of ^ as 
the parameters e vary; in particular, we want to prove the continuity 
of the solutions with respect to e in a suitable Banach space. 

To apply the results of the previous section we choose X = L^ and 
Y — H^. Before proving the main result, we collect some known facts 
about Sobolev spaces (in one space dimension). 

Lemma 1. (0) The map A'^ = (1—9^)2 is an isometric isomorphism 
of H^,s > into L^. The inverse of A^ is A~*. 

(i) Sobolev embedding (particular case): if u £ H'',s > -^ + n, then 
u is n-times differ entiable and there exists a constant c such that 

\\d>\\L^ < c Mj,. . (10) 

(ii) Algebra property: ifu,v e iJ'*, s > -^j then there exists a constant 
c(s) such that 

\\uv\\h.<c{s)\\u\\h.,\\v\\h. . (11) 

(Hi) Schauder estimate: i/a : M — > R is a smooth function and s > i 
then there exists a constant c{s^a, ||u||^a , ||w||^a) such that 

\\aiu)-a{v)\\^^<cis,a,\\u\\H^,\\v\\H^,)\\u-v\\^^ . (12) 

(iv) Commutator estimates: let u E H'',s > |. Then the operator 
Tu = {A^udx — udxA'')A^'^ is bounded on L^ and there exists a 
constant c{s) such that 

\\Tuh^<c{s)\\u\\^. . (13) 

Proof. For {iv) see [T3], Lemma A. 2. For all other statements, see [H] 
Appendix A. D 

We prove the following 



10 



Theorem 6. Let W he a hall of radius R in H'^ , such that u(0) G W . 
There is a time T such that for any e G K" the Cauchy prohlem (0) 
has a unique solution 

u{t-e) e C([o,r],p^)nci([o,T],i?"-(2»+i)). 

The map e ^ u(t;e) from W to C{[0,T],W) H C\[0,T], H''-^'^''+^'>) 
is continuous. 

Proof. As a first step we prove that u{t; e) exists and is continuous as 
a map from R" to C([0, T],i?-') n C^{[0,T],L'^). To prove this, it is 
sufficient to show that aii conditions of Theorems |4] and \E\ are satisfied 
uniformly in M", assuming X = L'^{R) and Y = H'{R). 

(X) It is trivially satisfied since X, Y are Hilbert spaces. The required 
isometry can be chosen to be A''. 

(Al) Following Example [2] above, we see that A{u; e) G G{X, 1, 13{R)), 
with 

P{R) = - sup sup |i9:Ea(u)| < -R sup a'{x) . 

2 new xeR 2 |a:|<i?, 

(A2) Since A''' commutes with the derivative operator, we have 

i?Ke) := (AM(e) - A(e)A^)A-^ = T,(„), 

where T is defined as in Lemma [T] (iv). Due to the commutator 
estimate and the Schauder estimate, we have that 

||-B(u;£)|i^.(R) <c{s,a,R)R^Xi{R). 

(A3) A{u;e) is a continuous operator from TP to H\ ioi < I < 
s — {2n + 1). Indeed, X^iLi EjC^k*^^ i'' continuous from H^ to H^ 
and a{u)dx is continuous from iJ'* to H"^^ due to the Schauder 
estimate and the algebra property of Sobolev spaces. A simple 
computation shows that 

\\{A{u; e) - A{v; e))\\jjs ^2 < c(s) sup \a'{x)\ \\u ~ v\\^2 , 

\x\<B. 

for some constant c{s) depending on s only. We choose 
/ii(i?) = c{s) sup |a'(a;)| . 

\x\<B. 

(A4) The same reasoning as in (A2) above, shows that property A(4) 
is satisfied with fJ.2{R) = c{s, a, R) for some constant c{s., a, R). 

(CI) The operator A{u;e) depends continuously (in norm) on e. In- 
deed, if e' ~ {e[, . . . , e'„), we have 

n 

\\{A{u; e) - A{u; e'))!!^^,^. < Y. 1^' - <\ ■ 
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(C2) The operator B{u; e) does not depend on e. 
(C3) The initial data do not depend on e. 

We have thus proved that there exists a T > such that u{t; e) exists 
and is continuous as a map from K" to C([0,T],if^) n C^{[0,T],L^), 
we want to prove that it is continuous also as a map from R" to 

C([0,r],iJ'*)nCi([0,r],iJ^-(2«+i)). Now the time derivative of m sat- 
isfies 

ut{t;e) = A{u{t); e)u{t; e), 

and A{u{t); e), for fixed t, e is a continuous operator from iJ" to i7'*-(2"+i) 
To complete the proof, it is enough to prove that the map 



A{,):W X 
is continuous. Indeed, 

\\A{v,e')y-A{u,e 



Jifs-(2n + l) 



< \\a{v)da:y - a{u)da:y\\j^, -(2^+1) + 



Y.i^'-e.)d, 



2i+l 



jys_(2„+l) 



By Cauchy-Schwarz and the Schauder estimate we have 

\\a{v)d^^y ~ a{u)d^^y\\fjs-(2n.+i) < c{s,a,R) \\u-v\\-^, WvWh- ■ 
Moreover ||Er=i(e^ " e^^^^+'^L.-c^^+i, < Er=i k^ -e.l IMh^- D 

Remark 2. Following U^, it is possible to prove a slightly stronger 
version of this theorem. Indeed, we can prove that the Cauchy problem 
(PI is uniformly locally well-posed in H^ with s > |. However, this 
stronger result is unnecessary for the purpose of studying the deriva- 
tives (with respect to £i) of the solution of the Cauchy problem ^ . 

Theorem [B] establishes that, fixed the nonlinearity a{u) and the ini- 
tial datum (/3 (E H^ , there exists a time T > such that the Cauchy 
problem ^ is locally well-posed in the time interval [0,r], continu- 
ously with respect to e 6 K" . In particular the life-span of the solution 
can be chosen independently on e. 

The natural problem is to find the supremum of all the positive 
times such that the Cauchy problem is locally well-posed and continu- 
ous with respect to e G C/ for some open U C K". We denote this time 
Tu. Suppose that for some a{u) and some Lp, there exists a J7 C M", 
such that the Cauchy problem is globally well-posed for all e E U. In 
this case, it follows that Tjj = oo. For example, if n = 1, ei ^^ 0, and 
linic^oo^^^ ~ 0, then the Cauchy problem is globally well-posed in 
iJ". for s > I [T3|. On the other hand, if for instance a{u) ~ u^, some 
solutions do blow-up at a finite time J19j . 

Since the general pattern is unknown, here we analyze the KdV 
Cauchy problem 



Ut 



eur, 



u\t=Q = ipe H', s > 3. 



(14) 
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If £ = 0, the solution of above Cauchy problem develops a gradient 
catastrophe singularity at a finite time t = tc > 0. Here tc coincides 
with the supremum of the positive time t for which the solution of the 
Cauchy problem can be continued. Conversely, if e ^ the Cauchy 
problem is globally well-posed |14) . 

We now show that the solution of the Cauchy problem of KdV is 
continuous with respect to e in any time interval [0, T], with T strictly 
smaller than the critical time tc- This is a simple corollary of Theorem 

m 

Theorem 7. Let T he any positive time smaller than the critical time 
T <tc and letu{t;e) G C([0,T], iJ'*) nC'^{[0,T],H''-^) be the unique 
solution of the KdV Cauchy problem (J14p . Then e — >■ u(t\e) is a 
continuous map fromW to C{[0,T], H") Ci C^HO^T], H"-^). 

Proof. The proof follows from Theorem [5] and a standard continuation 
argument. D 

Corollary 1. Let s > 3, T be any positive time smaller than the 
critical time T < tc and let u(t; e) be the unique solution of the KdV 
Cauchy problem (|14p . Then 

lim^^o \\u{t;e) — v {t)\\ , uniformly in t £ [0,r]. (15) 

Here, w" is the unique solution of the Cauchy problem for the Hopf 
equation 

v',=v%l v\^,=^. 



4 The ^-expansion of KdV solutions 

This section is devoted to study the differentiability of solutions of 
equation ^ with respect to e. For simplicity, we consider in detail the 
case of KdV only; as explained in Theorem [S] below, there is no trouble 
in extending the results to the whole class ^. 

We show that if s > ZN + 3, then the solution of the Cauchy 
problem of KdV P^ is A''-times differentiable with respect to e in any 
time interval [0,T], with T strictly smaller than the critical time tc- 
Before our main theorem, state a technical lemma. 

Lemma 2. Let / : M ^ C([0,r],i75) and g -.M.^ C([0,T],i7i) be 
two functions, and take X ~ L^ , Y ^ H'^ , with s > 3. Then, the 
family of linear operators 

^(£i, £2, £3) = /(£i)5^ + 5(£2) + £35^ (16) 

satisfy the conditions (i) through (iv) of Theorem Moreover, if 
{61,62, s^^) — > (£i,£2,£3) is any converging sequence, then the sequence 
of operators A" = A{6i,62,6^) satisfy conditions (vi) and (vii) of The- 
orem\^ 
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Proof. The operator A{ei, 62,83) was considered in Examples [5] and 
[3] above, where it was shown that it belongs to G{X, l,f3 + f3') with 
/3 = i sup3.gR \fx\ and /3' = sup^jg^ \g\. The verification of conditions 
(i) through (vii) follows the same steps as the proof of Theorem [HI D 

Theorem 8. Let N = [s/3 — Ij and [0, T] be a time interval such that 
<T <tc. Let 

u-.M.^ c([o,r],iJ") n cn[o,T],i?''-3) 

be the map that associates io e G M the unique solution of the KdV 
Cauchy problem |J^p . Then, there exist and are continuous the maps 

u^^) : R ^ c([o,r], ij'-^*'') n ci([o,r], ij'*-3(fc+i))^ 

for k = I, . . . , N, defined as 

Fixed e G M and /c e N, 1 < fc < N , then the function u^'^' (e) satisfies 
the following Cauchy problem 



k-l 






(17) 



w(")|t=o(£)=0, (18) 

where the linear operator A^ is defined as 

A, = uie)d^ + u^is) + edl (19) 

and we use the convention u^'^'{e) = u{e). 

Note that p7p is a linear non-homogeneous differential equation. 

Proof. If iV = 0, then the first part of the theorem follows from Theo- 
rem [3 while the second part is empty. Assume iV > 1, we now prove 
differentiability. To this aim we introduce the difference quotient 

^W(e,M= "^'^^^'"^'\ 
and a simple computation shows that u'^'(e, h) satisfies the equation 

ui'\e, h) = A(£, h)u^^\e, h) + u,,,{e + h), u^^\e, h)\t=Q - , 

where A(e,/i) — u{e)dx + Ux{e + h) + ed^. In the limit /i — > 0, 
the above equation converges to (|17p . with fc = 1. This is a linear 
non- homogeneous equation, with forcing term Uxxx{s). From The- 
orem [71 we have that Uxxxis) is a continuous function from M to 

c([o,r],/j^-3). 

Hence, we can prove the convergence of lim/i_>.o w^(£, ^) using the 
perturbation Theorem[31 provided that: i) we look for solutions of PT| 
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lying in the same space of the forcing term, namely H^~^, and ii) the 
forcing term belongs to D{A^). The latter condition holds since s > 6 
by hypothesis. 

Due to Theorems [2] and [3] and to Lemma [21 we conclude that 

u'^^He) := limu(i)(e,/i) 

solves (dH) and maps M continuously to C([0, T],H'-^)nC^{[0, T],L'^). 
Moreover, we have that the function w*^^) maps M continuously into 

C{[0,T],H'-'^)nC\[0,T],H'~^). 

Indeed, u). ' (e) equals A^u^^^^ + Uxxxi^), and the operator A^ maps 
any continuous function R ^> C{[0, T], H''~^) to a continuous function 
M — > C([0,T],_ff'*~^). The last statement can be proved in a similar 
way as in the proof of Theorem [6] 

We continue the proof by induction on the order of the derivative. 
Suppose the thesis is valid for i = 1, . . . , fc < iV. As before, we define 
the difference quotient 

that satisfies the non-homogeneous linear equation 

ui''+'\e, h) = A{e, h)u^'+'\e, h) + 41 (e + /i) + f+\e, h) 

+u^'\s)ui^\e)+u^:\s)u^^\e), (20) 



X: Q)u(^)(£)4^-^)(£) + fc4t^'(e) J . (21) 

The non-homogeneous term of equation ((20)) belongs to C([0, T],i7'*^^^'^+^^), 
and it continuously depends on e. In the limit /i — ;► 0, the quantity 
f'^^^{e, h) converges in C([0, T], i/''*^'^^'^'*'^'), continuously with respect 
to e, to 

Hence, the same reasoning as in the case of u'-^-'(£) shows that 
u'^'^+^Me) -.^ limu('=)(e,/i) 

h— s-O 
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solves the equation 

«f +^)(£) = A(£)u('^-+i)(£) + 45je) + .ft+^e) 

+iL^'\s)ui'\s)+ui'HsW'^\s), (22) 

and it maps M continuously to C([0, T], iT^-3('=+i))nCH[0, T], iJ^-3('=+2)). 
It is a simple computation to show that ([^^ coincides with p7|) . D 

Theorem [8] shows that if the initial datum of the KdV equation lies 
in the Sobolev space iJ", then the solution of the Cauchy problem is 
A^— times differentiable with respect to e, for N = [.s/3 — IJ . Con- 
sequently, if the initial datum lies in all the Sobolev space - e.g. it 
belongs to the Schwartz class - then the solution of the Cauchy prob- 
lem is smooth with respect to e. In particular, the solution admits an 
asymptotic expansion in power series of e. More precisely, we have the 
following corollary of Theorem [51 

Corollary 2. Let ip G H°° ~ ns>oH'' , T > be any positive time 
smaller than the critical time T < tc and u : M — > C([0, T], H°°) be the 
solution of the Cauchy problem {14^ . Then 

(i) u : R ^ C([0,r],if°°) is a smooth function (ofe). 

(ii) In e = 0, u admits an asymptotic expansion in power series of e: 



uis)^J2^'^'- (23) 



fc=0 

where v'^ — u{0) is the solution of the Cauchy problem 

v^ = v%l (24) 

«°|t=o=^, (25) 

for the Hopf equation, and v^ = u^^' (0) is the solution of the k-th 
linear Cauchy problem il'/]) when e = 0, that is: 

w1t=o = 0, fc>l. (27) 

Note that a similar Theorem was proven for the defocusing Non- 
linear Schrodinger equation |llj . 

Below, we state the analogue of Theorem [5] for the general equation 
© and we give a sketch the proof. The details of the full proof, which 
is rather long, will be given elsewhere. 

Theorem 9. Let U be an open subset o/M", [O,^] be a time interval 
such that the Cauchy problem is locally well-posed and continuous with 
respect to e & U and let 

u : M ^ C([0, TIH') n C^{[Q, r], iJ^-^^n+i)) 
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be the map that associates to e E U the unique solution of the Cauchy 
problem for (^. Moreover, let K = YJ-=i N^{2i+\) and N ^ X^ILi ^i- 
If s — K > 2n + 1, then the partial derivative 

d^u{e) 



exists m C([0,T],i7'^--^) nCi([0,T],if"~-f^-(2n+i)-) ^^g -^ continuous 
with respect to e E U . 

Proof. The Theorem can be proven along the very same hnes of the 
proof of the analogue Theorem [5] for KdV. More precisely, it is possible 
to prove existence and continuity of the partial derivative, by showing 
that it satisfies a linear non-homogeneous equation. Note that the 
condition s — K > 2n + 1 implies that the forcing term belongs to the 
domain of 9^"+^. D 



5 Hamiltonian perturbation of quasilinear 
conservation laws 

We now consider the results of the previous sections in the setting of 
the general construction, proposed by Dubrovin and Zhang [51 [^I10|. 
of Hamiltonian regularization of the quasilinear conservation law: 

ut = a{u) u^, u\t=o = f, (28) 

where a and the initial value (p are assumed to be smooth functions, 
and if is either periodic or rapidly decreasing at infinity. We discuss 
the aspects of the Dubrovin-Zhang construction which are more related 
with the present paper; in the next section we will show how the results 
obtained in Section [3] and 3] provide a rigorous justification to this 
method for a particular class of equations of type ^ . 

Let us consider equation ([^5)) . It is well known that this equation 
can formally be written as a Hamiltonian system 

SjTT 

with Hamiltonian 

H— h{u{x))dx, h"{u)~a{u), 



and where the Euler-Lagrange operator is defined as 
6H _ V-. ^fc rf'' dh 

for any local functional H ~ J h{u, u^, Uxx, • ■ • ) dx. The above Poisson 
bracket is given by 



err rjx [ ^^^ P, ^^^ A 

{Hi,H2}= / T-r-^9xj—-:dx, 
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for any pair of functionals Hi ^ J hi{u, u^, Uxx, ■ ■ ■) dx, i ~ 1,2. 

FoUowing Dubrovin, by Hamiltonian regularization (or perturba- 
tion) of the quasilinear conservation law (|28p we mean an expression 



Ut = \u,h\ = dx^—r-r, u\t=o = ip, (30) 

L J duix) 



SH 

' Su{x) 

where the Hamiltonian is given by a formal series 

H = H + Y,Hke\ Hk^ f hkiu;ux,...,ui^'^)dx, fc > 1, 

k>l •' 

for some hk, which are assumed to be differential polynomials in the 
derivatives. In addition, the solutions of equation pop are sought to 
be of the form 

M(x,i)=^w*(a;,i)e2i, (31) 

i>a 

with coefficients u* smooth functions of x and t. Within this setting, 
all identities are understood in the sense of formal power series in e 
- they are assumed to hold identically at every order in e. Therefore, 
the perturbed Hamiltonian and the solutions are not required to be 
convergent (neither asymptotic) series. Note, however, that the initial 
values of the perturbed and the unperturbed equations are assumed 
to be the same. In particular, the function tp is independent of e, and 
from the expansion ([31]) we deduce the identities 



v°{x, 0) = (p{x), v'{x, 0) = 0, Vi > 1. (32) 

A primary task in the Dubrovin-Zhang approach is the classifica- 
tion of Hamiltonian perturbations, which is performed by considering 
equations pop modulo quasi-Miura transformation. These are trans- 
formations of the form: 

u^^« = ^e''^Ffe(u;M„...,uW), (33) 

fe>0 

where the functions F^ are rational with respect to the derivatives. A 
partial result is given by the following 

Theorem 10. J^ Any Hamiltonian perturbation of equation (j28p of 
order e* can be reduced by a Miura-type transformation to an equation 
of the form (|30p , with Hamiltonian 



H= h{u;Ux,Uxx,s)dx, h = a, (34) 

h^h-'-^ch'"ul+e^\^(^ph'" + ^c^h(^^yi,, 

- (^ /.(^) + ^ h(^^ + 1^;^(«) + P-hi^^ + ^/^(^) - . /.'") ut 

for arbitrary functions c{u), p{u), s{u). 
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Let us now consider in more detail the solutions of the perturbed 
equation pop , which - after Theorem [10] - we will consider together 
with a Hamiltonian of the form ([M|) . By expanding both sides of 
equation pOp according to the Ansatz pip , in first approximation one 
obtains 

v^ = aiv')vl (35) 

which says that v'^ must be a solution of the unperturbed equation (pS)) . 
Accordingly, from the higher order coefficients one obtains an infinite 
set of linear non-homogeneous equations (or transport equations) for 
the coefficients v^{x,t), which can be solved recursively starting from 
the solution of ((35)) . For instance, the equation for v^ turns out to be 

vl=dJav^^ca'vl^^'^{ca" + c'a'){vlf\ v\^o ^ (36) 



Remark 3. In the class of initial data considered in the present paper, 
any solution of equation psp develops a singularity at a finite time 
t ~ tc > 0, known as time of gradient catastrophe. An important 
aspect of Dubrovin's theory is concerned with the study of the solution 
of pop in a neighborhood of the critical time tc, in order to show how 
the singularity is regularized by the dispersive perturbations \^. 

This, however, is out of the scope of our present method of investi- 
gation. In what follows we study the series pip for any time interval 
[0,T] strictly smaller than the critical time: T < tc- 

Example 4. The KdV equation 

Ut=U U:c + e^ Uxxx 

is obtained from pOp . p4p by choosing h{u) = g"^; c{u) = 1, and 
p{u) = s{u) = 0. In this case, a simple computation shows that the 
Cauchy problems for the transport equations are given by 

^t=il{%'<~'+kvl-l k>0, 

One can - in principle - solve these equations recursively. Note that 
these Cauchy problems coincide with p4p . (j26p . 

From the above discussion it follows that - at least in principle - 
all coefficients of the expansion (PT|) can be obtained once a solution w° 
of equation (|55)) is known. The method followed in |F1 IH] (see also the 
older result [Tj) to find solutions of the perturbed equation (pO)) is to 
construct a quasi-Miura transformation, relating the solution v^{x, t) of 
the unperturbed equation ([55)) to the the solution u{x, t) of perturbed 
equation ([50)) . The required transformation has been suggested in [S] 
to be of the form: 



^0 ^ ^ ^ „o 



{v\x),K]+e^{{v\x),K],K] + ... (37) 
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where the functional A', up to order 4 in e, is given by 



K 







This in particular implies 



v' = \d^[c{v'f^+c'{v')vl), (39) 

and a direct substitution shows that the above function satisfies equa- 
tion ((55|) . provided v^ satisfies ([55]) . Note, however, that the function 
((39|) is bounded only for monotone solutions of equation ((35l) . Further- 
more, (|39p does not satisfy the required initial condition 



6 Solutions of the transport equations 

The classification problem of Hamiltonian perturbations, together with 
the quasi-Miura transformation discussed above, are main ingredients 
of the Dubrovin-Zhang constructions before the critical time, which is 
the time-span we are interest in the present paper. As already noticed, 
this approach is mainly based on identities of formal power series, and 
indeed, if for a given Hamiltonian the corresponding equation pop is 
just a formal series, then the only way to construct a (formal) solution 
seems to be through the use of a formal series, for instance like ([5T|) . 

However, if the Hamiltonian perturbation (|30p is well-defined - for 
example if the Hamiltonian H is given by (j34p - then the resulting 
equation may happen to be locally well-posed in some function space, 
say 7J^, for s big enough. At the same time, the solution may be 
A^— differentiable with respect to e and the formal series (PU)) may just 
be - up to order A^ - the Taylor expansion of the actual solution. 

Using the results of Section [3] and 2] we can show that this is true 
if the Hamiltonian perturbation ([50]) coincides with a generalised KdV 
equation ^. This is the case for the general equation Q, provided 
n < 2. Indeed, the following Lemma holds: 

Lemma 3. The equation 

ut = a{u) u^ +e^ a u^^^ + e^ (i u^xxxx, (40) 

which is obtained by ^ in the case n — 2 and ei — a e"^ , 62 = (3 s'^ , co- 
incides with the Hamiltonian perturbation llclU\) with Hamiltonian p4p . 
provided the coefficients are chosen in the following way: 

a f3 3a'^ a" 

2 f2 (a")^ 7 a" a'" , 1 a"" ^^ /3 f [a")'^ 



5 (a')5 20 (a')4 24 (a') V 12 V (a')^ («')^ 
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Theorem 11. Let ip G U^ ., s > 11, U a neighborhood of e ~ 0, u{e) be 
the solution of the Cauchy problem (|40p with initial data ip and [0, T] 
be a time interval such that the Cauchy problem is locally well-posed 
for any £ G [/ and continuous with respect to e G U. Moreover, let 
s > 5 + 6A^. Then u{e) has a Taylor expansion in s = 0, up to order 
AN , of the form 

2N 
fc=0 

where 

v^ G C([0,r],iJ'*-3fc)^ ^(g) g C{[Q,T],H'~^^). 

Here v^ = u(0) is the solution of the Cauchy problem for unperturbed 
equation (j35p . while v^ satisfies the first transport equation (j36p with 
parameters (j4ip . v^ ,k > 1 is the solution of the higher transport equa- 
tions. The remainder term r is o(e*^). 

Moreover, if a = then the same is true, under weaker hypothesis. 
Namely, ifs>5 + 5N, then v'' ^ for k odd, v"^' G C([0, T], iJ'*-^'), 
andr -.R-^ C{[0,T, H''^^]) is continuous, with r = o{e'^^). 

Proof. The proof follows from Theorem [HI D 

Remark 4. Due to Theorem Tg > for any Q C M. //, for 
(e, a, /3) ^ (0,0,0) the Cauchy problem is globally well-posed }16^ . then 
the time Tq in the above theorem can be chosen to be any positive time 
smaller than the critical time tc of the unperturbed equation i35]) . This 
is the case, for instance, if a{u) = u. 

In the following theorem we provide a correction to the formula 
([55)1 for the first coefficient v^ . Om- formula turns out to be valid 
for solutions of (j35p which are not monotone, and satisfies the correct 
initial value. 

Theorem 12. Consider the Cauchy problem for equation (|35p with 
initial datum i;*'(0) = </? G H^ , s > 3, denote by tc be the associated 
critical time, and let v'^{x,t), with (x,t) G M x [0,tc), be its unique 
classical solution. Then, the solution of the linear transport equation 
(|36p . with v^ given above and initial datum v^{{)) ^ is 

where the family of functionals Kt, with t G M, is defined by 



KM--^-^ 



- I c(u)uj,\og{l + ta {u)ux) dx, 
for every u G H" , with ||m||^s small enough. 
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Proof. The explicit form of the function v^ stated in the theorem is 
given by 

,_t d ({ca'nvlf + 2ca'vl,+tc{a'fvlvl,+tc'{a')\vir 



V 



2dx\ (l + ia'w0)2 

(43) 
where the functions a, c, and the corresponding derivatives are eval- 
uated at u = v'^{x,t). A direct calculation shows that this function 
satisfies equation (P5|) with the correct initial value. D 

Remark 5. Note that the above theorem holds for any choice of the 
functions a and c. This fact suggests that similar results of the one 
obtained in Theorem \ll\ remain true for a generic Hamiltonian pertur- 
bation of the quasilinear conservation law (|28p . 

Remark 6. Formula (|43p has been obtained making use of the so called 
'string equation', introduced in the setting of Hamiltonian perturbations 
of nonlinear PDEs by Dubrovin JBj. Although heuristic, the use of the 
string equation turns out to be a very powerful method for describing 
solutions of the perturbations both before the critical time and in a 
neighborhood of it. 
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